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Abstract 

A review of the appearence of integrable structures in the matrix model desrip- 
tion of 2d-gravity is presented. Most of the ideas are demonstrated at the techni- 
cally simple but ideologically important examples. Matrix models are considered 
as a sort of "effective" description of continuum 2d field theory formulation. The 
main physical role in such description is played by the Virasoro-VF conditions, which 
can be interpreted as certain unitarity or factorization constraints. Both discrete 
and continuum (Generalized Kontsevich) models are formulated as the solutions to 
those discrete (continuous) Virasoro-VF constraints. Their integrability properties 
are proved, using mostly the determinant technique highly related to the representa- 
tion in terms of free fields. The paper also contains some new observations connected 
to formulation of more general than GKM solutions and deeper understanding of 
their relation to 2d gravity. 
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1 Introduction 



The most geometrical way to derive the 2d gravity or string theory partition function (as 
well as the generating function for the correlators) is given by the Polyakov path integral 
which for the partition function reads 

F(X)= £ X p F p 

p— genus 

F p = £ Dgexpj J RA-'R (1) 

Matrix models historically appeared |T| when one considers the discretization of ([I]) 
which in a sense turns to be exact at least for the simplest cases of "empty" string theories, 
when the target space is low-dimensional (or in the limit of pure gravity does not exist at 
all). A well-known example of such theories is given by minimal (p, q) models coupled to 
2d gravity. 

The difficulties with the continuum formulation (|1|) as usual arise from the fact that 
it possesses additional complicated information (like the "Verma- module" structure of 
the underlying 2c? CFT) which might not be essential for description of final "effective" 
theory. We have exactly here the case of a gauge-theory, when after coupling to gravity 
all conformal descendants are "gauged away" and one might search for a sort of "frame" 
description, which luckily appears in the form of matrix models. 

Matrix models are usually defined by the integrals like 

Z N = J DM NxN exp[-TrV(M)] (2) 
which in the continuum limit, requiring in part iV — > oo 

logZ - F (3) 

can give the whole nonperturbative solution to ([!]) @. 

However, below we will mostly advocate the point of view different from the original 
definition (|2]). Indeed, the main difference between the continuum ([!]) and the matrix 
model formulation is that the first one requires some sort of the unitarity or factorization 



model formulation (0) these relations usually appears automatically. Moreover, for known 
solutions they usually appear in the well-known form of the Virasoro-W 7 [] constraints, 
which in a sense may be considered as a definition of matrix model . 

Below, we will usually start the description of various matrix model and matrix model 
-like theories as being particular solutions to these Virasoro-P^ recursion relations. It 
turns out that these relations lead to the integrability property of matrix models, namely 
the solutions to these constrains turn to be r-functions of the hierarchies of well-known 
integrable equations || £|, ||]. 

In terms of the partition function ([[J) , (|3]) 

F(T) = log r(T) (4) 

where T = {T^} is the set of times or the coupling constants in the theory of 2d gravity 
The formula @ or the appearence of the integrability is exactly what allows one to 
make more progress in studying "frame" formulation (0) instead of the original one ([I]). 

In what follows, we will first consider an example of discrete matrix models (finite N 
in (H)) as being solution to the most simple discrete Virasoro constraints and then pass to 
the continuum case. Both kinds of solutions to the recursion relations are r-functions of 
well-known hierarchies and both particular solutions have representation in integral form. 
We will also stress the moment that the solutions to discrete constraints correspond to 
the discretization of the world-sheet of string, while, as it will be seen below the matrix 
solutions to the continuum relations have rather interpretation of target-space theory or 
of an effective string field theory. 

1 This should have an interpretation as a sort of world-sheet - target-space duality between the Virasoro 
symmetry as gauge symmetry for 2d gravity and Virasoro relations in the target-space (or better in the 
space of coupling constants). 



2 Discrete hermitean 1-matrix model as a solution to 
the discrete Virasoro constraints 



In this section we are going to consider the first (and simplest) example, demonstrating 
the above ideas, namely the solution to the discrete Virasoro constraints ||: 

L n Z[t] = 0, n > -1 

oo 

L n = J2kt k d/dt k+n + ]T d 2 /dt a dt b (5) 

k=0 a+b=n 

with an additional requirement (concerning £ -variable) 

dZ N /8t = -NZ N 

(later on N will be identified with the size of matrices in the formulas like (0)). The key 
idea how to solve the constraints ([5]) appears after one notices that the Virasoro generators 
(IH) actually have the well-known form of the Virasoro operators in the theory of one free 
scalar field 0. If we try to look for such solution in terms of holomorphic components of 
the scalar field 

4>(z) =q + p\ogz + ^2^z~ k 

k^O fc 

[J n ,Jm] = n5 n+mfi , [q,p] = 1 (6) 
the procedure is as follows. First, define the vacuum states 

J k \0) =0, (N\J_ k = 0, k>0 

p\0) =0, (N\p = N(N\ (7) 
"Half of the stress-tensorQ components 

T(z) = \[d<P{z)f = Y J T n z~ n -\ T n = Wj-kJk+n + \Y, J ^ ( 8 ) 

a.b>0 



2 In this sense we have c = 1 for (g) though it is not too sensible to speak of the central charge, having 
only the half (n > — 1) of the Virasoro algebra, i.e. 

[L n , L m ] = {n- m)L n+m , n, m > -1 

always without the central element. 

3 For the sake of brevity, we omit the sign of normal ordering in the evident places, say, in the expression 



obviously vanish the SX(2)-invariant vacuum 

r„|0) = 0, n>-l (9) 
Second, we define the Hamiltonian by 

H ® =4;E«= / V(z)j(z) 

V(z) = £t^, j(z) = ^d<t>(z). (10) 

fc>0 v l 

Now one can easily construct a "conformal field theory" solution to (Bj) in two steps. The 
basic "transformation" 

L n (N\e H V... = (N\e H ®T n ... (11) 
can be checked explicitly. As an immediate consequence, any correlator of the form 

(N\e H ®G\0) (12) 

(N counts the number of zero modes, "included" in G - that is the role of the size of 
matrix in @) gives a solution to fl5|) provided by 

[T„,G]=0, n>-\ (13) 

The conformal solution to ([13]) (and therefore to (|5|)) immediately comes from the basic 
properties of 2d conformal algebra. Indeed, any solution to 

[T(z),G]=0 (14) 

is a solution to (|i~3D , and it is well-known that the solution to ( |i~4l) is (by definition of the 
chiral algebra) a function of screening charges in the free scalar field theory given by 

Q ± =<fj ± =<fe ± ^f. (15) 



With a selection rule on zero mode it gives 

G = exp Q + - ^Q N + (16) 

(Of course, the general case might be G ~ Q_ but the special prescription for 

integration contours, proposed in 0, implies that the dependence of M can be irrelevant 
and one can just put M = 0.) In this case the solution 



after computation of the free theory correlator, analytic continuation of the integration 
contour gives the result 

r N 

Z 2 ,n = (Nir'j n^exp(-X)^) = 

i=l 

= (NNo\ U(N))- 1 J DM exp (- J2 tkM k ) (18) 

N 

A N =Y[(Zi-Zj) 

i<j 

in the form of multiple integral over the "spectral parameter" or, in this particular case, 
in the form of the integral over Hermitian matrices of the type of eq.(^). 

This point of view actually could be considered as a constructive one. Namely, instead 
of considering a special direct multi-matrix generalization of flTSP one can use powerful 
tools of 2d conformal field theories, where it is well known how to generalize almost all the 
steps of above construction: first, instead of looking for a solution to Virasoro constraints 
one can impose extended Virasoro or jy-constraints on the partition function. In such 
case one would get Hamiltonians in terms of mn/tz-scalar field theory, and the second 
step is generalized directly using screening charges for VT-algebras. The general scheme 
of solving discrete jy-constraints looks as follows |7]]: 

(i) Consider Hamiltonian as a linear combination of the Cartan currents of a level one 
Kac- Moody algebra Q 

...,^ e ))= £ 4 A W,, (19) 

A,fc>0 

where {/J-i} are basis vectors in Cartan hyperplane, which for SL(p) case are chosen to 
satisfy 

1 p 
P j= i 

(ii) The action of differential operators with respect to times {t^} can be now 
defined from the relation 

wt\N\e Hi - {t}) ... = (iV|e" ({ ' }) W; a) ■ ■ ■ , a = 2, . . . ,p- i > 1 - a, (20) 

where 




are spin-a W-generators of W p -algebra written in terms of rank ^-component scalar fields 

(iii) The conformal solution to the discrete ^-constraints arises in the form [0 

^ C iV [{*}] = (N\e H ^G{Q {a) }\0) (22) 

where G is again an exponential function of screenings of level one Kac-Moody algebra 
(see and references therein) 

QW = ^jW = ^ e ^ (23) 



{a} being roots of finite-dimensional simply laced Lie algebra Q. The correlator ( p^) is 
still a free-field correlator and the computation gives it again in a multiple integral form 



^ P c iv M [«] ~ /n 



i=l \ A,fc>0 



X 



x nniK^-^r (24) 

(a,/3) i=lj=l 

The only difference with the one-matrix case (|T8| ) is that the expressions have rather 
complicated representation in terms of multi-matrix integrals, the following objects will 
necessarily appear 

ft 8(4°° - z f ) )° L ' 3 = [det{M^ <g> I — I ® M^}] al3 , (25) 
i=i j=i 

However, this is still a model with a chain of matrices and with closest neibour interactions 
only (in the case of SL(p)). 

Actually, it can be shown that CMM, defined by (p2|) as a solution to the PF-constraints 
has a very rich integrable structure and possesses a natural continuum limit |7], [TO]. To 



pay for these advantages one should accept a slightly less elegant matrix integral with the 
entries like (p5|) . 

The first non-trivial example is the p = 3 solution to H/3-algebra: an alternative to the 
conventional 2-matrix model. In this case one has six screening charges Q^ ±a ^ (i = 1,2,3) 
which commute with 



and 

W^(z) = j:(» x d ( f>(z))\ (27) 

A=l 

where /x A are vectors of one of the fundamental representations (3 or 3) of SL(3). 



The particular form of integral representation fl24j) depends on particular screening 
insertions to the correlator (p2|) . We will concentrate on the solutions which have no 
denominators. One of the reasons of such choice is that these solutions possess the most 
simple integrable structure, though the other ones can still be analyzed in the same 
manner. 

The simplest solutions which have no denominators correspond to specific correlators 

2ffi[{t}] = (N\e B «*»l[e*pQ ai \0) (28) 

i 

when we take a; to be "neighbour" (not simple!) roots: (c^ck,) = 1. In the case of SL(3) 
this corresponds, say, to insertions of only Q ai and Q a2 (again, we use the notations of 
U) and gives 

Z°™? N [t,t\ = Z MjN [t,t\ = ^(N,M\e H ^(Q^f(Q a2 ) M \0) = 
= T^rry / n^< d W«P ("E^te) + V{vi)]) A%(x)A 2 M (y) H( Xi - Vj ) (29) 
Possible generalizations of the above scheme can include the " supersymmetric matrix 



models" [11], where the authors looked for a solution to the system of equations L n Z = 
and G m Z = 0, with generators {L n , G m } forming the N — 1 superconformal algebra, 
which in our language is nothing but a trivial generalization of the one-field case which 
has to be substituted by a scalar superfield. Then the insertion of screenings of N — 1 
superconformal algebra immediately leads to the result of [pTTJI . From this point of view 
the real problem with supersymmetric generalization can arise only in the N = 2 case 
because of the lack of appropriate screening operators. 



2.1 Determinant representation and integrability of the solu- 
tions to Virasoro- W constraints 

In the simplest case of the Hermitean one-matrix model the N x N matrix integral flT8| ) 



= (Pi,Pj) = / Pi{m)Pj{m)e- v{m) dm = dye*'® (30) 



and equals 



N-l 



Z v I [ e^ {t > (31) 



i=0 



It follows from (j30| ) and the definition of orthogonal polynomials 



Pi{m) = ^dijin 3 



(32) 



an = 1 



that 



diag(e^ {t) ) = AHA T , (33) 
where A = \\a,ij\\, A T - transponed matrix, and H is so called matrix of moments 

Hij = J m i+J e- v(m) dm (34) 

Thus, 

Z N [t] = det[diag(e^ (i) )] = det AHA T = det H — r N [t] (35) 
Without going into all the details, which could be found in |12| , |i~3| we will only point 



out that (^) is a kind of a Riemann-Hilbert problem and the determinant formula ( |35"D 
is one of the basic definitions of r-function of Toda theory. 

The relation ([H]) and (|35|) is actually based only on the fact that in the theory of 
Toda-chain there exists a relation between the potentials of the Toda-chain equations and 
the r-function having the form of a difference operator: 

e w(t) = W!) (36) 

T i\t) 

Now to check that we have really got Toda chain hierarchy, let us check the first 
equation - flow in the direction This might be done after we introduce the Lax 



operator for the Toda chain (see |12[ for details), which in the basis of the orthogonal 
polynomials (|30"D acts by: 



mPiim) = P i+1 {m) - Pi(t)P(m) + R^P^m) (37) 

i.e. is determined by a trilin ear matrix, what follows from a sort of "unitarity condition" 
or just the properties of the basis (^). Now, from (]30f) , (|37D one can easily establish the 
relations among the "potentials" {(pi(t)} and the matrix elements {Ri(t) , Pi(t)} , first: 

(i\m\i - 1) = e w( *> = Ri(t)e^- l{t) (38) 

gives 

i?.(t) = e Vi(*)-w-i (39) 
Now, differentiating (|30|) for i = j one gets 

J-(^li) = e^ (t) ^ = / dme-^ m \-mPf + 2P,^) = Pi (t)e^ (40) 
ati OTi J oti 

where the second term in brackets dissappear again from the orthogonality condition and 
property (^). From ([E]) it follows that 

«W = ^ (41) 

is just the momentum for the </?j(t)-coordinate. 
Differentiating (|30|) with z > j, we obtain 



(42) 



= - y (fme-S^ 4 fmP;Pj + P,^ 
comparing which with fl3"T|) and using (RO) one gets 



9*i 

Now we are ready to differentiate 



dP 

' ' RiPi-i (43) 



dp 9P,-4-i cto* „ 



Multiplying ([0]) by Pi integrating and using (|30| ) and ( [43"D one finally gets 



= - J* (45) 

OTi 

or using (gTj), (gDJ) 

_|£ = gw+i-w _ e w-w-i ( 46 ) 

which is nothing but the first Toda-chain equation. In terms of the r-function ([46]) can 
be rewritten in the Hirota bilinear form: 



r N (t)^r N (t) - (^P) =T N+1 (t)T N - 1 (t) 
The r- function of p = 2 case (|35"D can be also written in the form 



(47) 



Z 2>N (t) = det N [d i+ i- 2 C(t)) = r N (t) (48) 

with 

d tn C(t)=dlC(t) (49) 
Eq.(|4*9|) means that C(t) just has an integral representation 

C(t) = J d(jL(z)exp y E, t '° zk > ( 5 °) 

where dfi is some measure; the Virasoro constraints ([5]) fix the concrete measure (dfi = dz) 
and the contour of integration in (|50D . The determinant form (|48| ) is an explicit manifes- 
tation of the fact that the partition function does satisfy the Hirota bilinear relations, the 



simplest one of which in this particular case takes the form ([47]). 

Now one can generalize ([48]) and ( f4~9|) 0. In the case of the p = 3 model (|22"D we have 
to introduce two functions instead of (pH|): 

C{t) = J dz exp[-V(z)], C(t) = J dz exp[-V(z)\ (51) 

where 

k>0 k>0 

and 



Again, we can get the determinant representation, now having the form (d = d/dti, 
d = d/dh) 

C dC ... d N - 1 C C BC ... d M - l C 

dC d 2 C ... d N C dC d 2 C ... B M C 



det 



QN+M-l(j qN+Mq q2N+M-2(J qN+M-1q qN+M(J q2N+M~2(J 

= T N,M(t, t) 



(53) 



which is exactly the double- Wronskian representation of a r-function [ 14[ . 

From representation (|5^) it is easy to derive the analogs of the Hirota relation (|47j ) 



d 2 , , js T~N+l,M~l(t,t)T N _ 1M+ i(t,t) 
j- log T N> M{t,t) = 



"Higher-times" Hirota relations have more complicated form. 



(54) 



2.2 Fermionic representation 

Now we shall proceed to the representation of the solutions to the Virasoro-H 7 constraints 
in terms of the fermionic correlation functions [] developed for generic integrable systems 
in 



15 



Indeed, the r-function of 2-component KP hierarchy defined by the fermionic correlator 



tP m {x, y) = (N, M\e H ^G\N + M, 0) 



with 



k>0 

JM( Z ) = J2 Jk ) z- k - 1 = :^{z)^*{z): 

5i 



^\z)^>{z') 



z — z' 



+ 



(55) 

(56) 

(57) 
(58) 



Indeed, the two types of technique we are using are practically equivalent: the symmetry (better 
antisymmetry) of determinants under permutations is what reflects the anticommuting nature of the 
fermions 



is equivalent to (|T7| ) for certain G when ( |55| ) depends only on the differences Xk — Uk- To 
prove this we have to make use of the free-fermion representation of SL(2)k=i Kac-Moody 
algebra: 

j = - ^ {2 V (2) *) = \{J {1) ~ J {2) ) 

J+ = ^(D* j_=^y 2) * (59) 

Now let us take G to be the following exponent of a quadratic form 

G = : exp (J ^ (2 V (1) *) : (60) 



The only term which contributes into the correlator (p5|) due to the charge conservation 
rule is: 

1 / r \ N 
G NM = G Nt - N 5 M ,-N =w/\j ^ (2 V (1) *J : 5m,-n (61) 

Now one can bosonize the fermions 



jW = «90 1; J (2) = <90 2 (62) 

and compute the correlator 

r%\x,y) =T$_ N (x,y) = ^(N,-N\exp ( 5>*4 1} + vrf } ) J (/ r^ty^:) " |0> = 



^^-iV|exp(/[X(,)J«(,) + y(,)^(,)]) (/ : exp^x - 2 ) :) |0> 
Introducing the linear combinations y/2<p = 0i — 2 , v20 — 0i + 02 we finally get 

(x,y) = ^(exp f-L j[X{z) + Y(z)]d${z)\ ) x 

x(N\exp(^=f[X(z)-Y(z)]d(f>(z)\ (| : exp V20 :) " 1 0) = t$ (x - y) (63) 

since the first correlator is in fact independent of x and y. Thus, we proved that the 
r-function (p5|) indeed depends only on the difference of two sets of times {xk — Vk}, and 
coincides with flTTD . 

The above simple example already contains all the basic features of at least all the 



situation. In other words, the diagonal £7(1) GL(2)-current J = \{J^ 1 ' + J^ 2 ') = ^d(f) 
decouples. This is an invariant statement which can be easily generalized to higher p 
cases. 

In the case of SL(p) we have to deal with the p-component hierarchy and instead of 
for generic r-f unction one has 

T jp\ x ) = (N\e H{x) G\0) (64) 
N ={N U ...,N P }, x = {x (1) ,...,x (p) } 

and now we have p sets of fermions {-0W*,-0W} i — 1, ...,p. The Hamiltonian is given by 
the Cartan currents of GL(p) 

H{t) =EE4 J) 4 i} (65) 

8=1 fc>0 

and the element of the Grassmannian in the particular case of CMM is given by an 
exponents of the other currents 

= jW) = ^W^0>, jfo> = ^(«)*^(7> i ^ j (66) 

i.e. 

G = nexp(g fe) )exp(g^)exp(g^>) (67) 

Q{ij) = j Q(V) = j J(«) j Q(ij> = j jfe> ; j ^ j 

Since (^) are the GL{p)i Kac-Moody currents, ( |67| ) play the role of screening operators in 
the theory. It deserves mentioning that they are exactly the SL(p) (not GL(p)) -screenings 
and thus the r-function (|64]) does not depend on {J2i=i x k^}- 
In the case of SL(3) this looks as follows. The screenings are 

Q(«) = ^J(«), (68) 

where {a} is the set of the six roots of SL(3). In terms of fermions or bosons the screening 
currents look like 



(55 



J(CX 2 ) 


= ^ {2) *^ (3) * 


= exp(0 2 + 


h 


j(a 8 ) 


= ^ (1 V (3) * = 


= exp(0 3 - <f 




j(-ao 


= ^ (1 ty (2) = 


exp(-0i - 


02 


j(-a 2 ) 


= V 3) = 


exp(-0 2 - 


03 


j(-a 3 ) 


= v {3 V (1) * = 


= exp(0! - q 


>s) 



(69) 



The particular r-function is now described in terms of the correlator 



T ®(x) = (N\e H{x) G\0) (70) 

with 

G-HexpQ^ (71) 
a 

The condition of Cartan neutrality is preserved by compensation of charges between the 
operator ( |7I|) and left vacuum (N\ in (|70|). It is obvious that in such case due to the 
condition of Cartan neutrality of the correlator (like in Wess-Zumino models) the mode 
J = d(f) = -t= Ya=i dtfti decouples from the correlator, and 



T ®{x) =(N\e H ^G\Q) 



(0|exp E^E4° |0>^ H( ^G|0> 



v fc>0 i=l 



(72) 



where the first correlator in the second row is trivially equal to unity. For the specific 
choice of the operator G in ([72] ) 

G = G lj2 = exp (7 J^) exp (| J< a '>) (73) 

we reproduce the formula (f29|). 

Finally, let us only stress two main ideas we have demonstrated above: first, we proved 
that the solutions to the discrete Virasoro-W 7 constraints can be rewritten from (free) 
bosonic to (free) fermionic language, which means automatically that they are solutions 
to integrable systems in the sense of [O . Second, in general case the solutions to discrete 
constraints are presented in the form of conformal multimatrix models, being particular 
solutions to mtzftzcomponent hierarchies. 



3 A solution to the continuum Virasoro -W constraints 



This section is devoted to the derivation of the solution to the continuous Virasoro and W- 
constraints. To be more precise, we shall investigate them as being the direct consequence 
of the " matrix" equations, which could be treated later on as the Ward identites for certain 
matrix integrals. 

This is however not exactly the same what we had before for the case of discrete 
constraints. The reason is that the continuum case differs from the discrete one roughly 
speaking by replacement of ordinary free scalar fields by the same fields but with antiperi- 
odic boundary conditions. This is the sense of so-called double scaling limit (see below 
for details) which can be done successfully for all conformal multimatrix models @, |i"0fl , 
discussed above. The exact solution in terms of conformal correlators is much more com- 
plicated for the antiperiodic fields f], thus instead here we are going to reformulate the 
problem. 

Fortunately, it turns out that the continuum Virasoro constraints can be summed up 
into certain matrix differential operators. Namely, for the Vy^-algebra these operators 
are related to the Laplacians (or better Casimirs) for corresponding algebras, having the 
form of 

+ ... (74) 



dAP 

where A is N x N hermitean matrix (for the SU(N)-case). These equations might be 
identified with the Ward identities for "continuum" matrix theories. 



So, we start with an equation of the type of (|74D operator vanish certain function, and 



prove that it is equivalent to the continuum Virasoro (or W) constraints. 

{Tr e(A)[V '(d/dA tr ) - A]}.F[A] = 0. (75) 

(|75|) can be certainly interpreted as a ward identity satisfied by a matrix integral which 
after the proper normalization gives a solution to continuum 2d gravity. The exact formula 
for corresponding partition function reads 

Z {N) [V\M] = C^ N \V\M]e TrV{M) ~ TrMV ' {M) ( DX e -TrV(X)+TrV'(M)X ( 7 Q) 



where the integral is taken over N x N "Hermitean" matrices, with the normalizing factor 
given by Gaussian integral 

CW\y\M\~ X = J DY e -TrU a [M,Y\ } 

U 2 = lim \Tr\V{M + eY) - V(M) - eYV'(M)} (77) 

and discuss only specific potentials, V(X) = const ■ X p+1 , giving rise (while substituted 
in (0)) to Q 

In the simplest example of p = 2 we have quadratic operator (or just Laplacian) and 
will prove the identity 

ifr(e^r - eA)F C " Z tr(eA- n - 2 ) (78) 

for 

( DX exp(-trX 3 /3 + trAX) = C[VX] exp(-trA 3 / 2 )Z^(T m ) 
J 3 

T m = —TrM~ m = —TrA~ m/2 , m - odd (79) 
m m 

C[VX] = det(VX®I + I®VX)-^ (80) 
and 



^ 2 >{A} = 



with 



1 ^ d 1 d 2 

2 fe>s„ + i ]0 k 9T k+2n 4: a +b=2n dT a dT b 

fc odd a,b>0 ; a, 6 odd 

T 2 1 9 

4 16 dT 2n+3 



While (|78 ) is valid for am/ size of the matrix A, in the limit of infinitely large A 
(N — > oo) we can insist that all the quantities 

6 The proof of the Virasoro constraints for generic potential is based on integrability and will be 
presented below. 



tr(eA- n - 2 ) (82) 

(e.g. trAP~ n ~ 2 for e = A p ) become algebraically independent, so that eq. ( |75|) implies 
that 

C n Z{T} = 0, n > -1 . (83) 

Note that ^"{A} in (|T9[) , which we have to differentiate in order to prove ([T8"D, depends 
only upon eigenvalues of the matrix A. Therefore, it is natural to consider eq.(78) at 



the diagonal point Ajj = 0, i ^ j . The only "non-diagonal" piece of (|78|) which survives 
at this point is proportional to 



d 2 X k 



dAijdAji 



= for i + J. (84) 

A mn =0, m^=n 1 J 



Eq . (]84|) is nothing but a familiar formula for the second order correction to the Hamil- 
tonian eigenvalues in ordinary quantum-mechanical perturbation theory. It can be easily 
derived from the variation of determinant formula: 



5log(det A) = tr^-SA - ^tr(-5A-5A) + ... . (85) 
For diagonal Aj,- = Aj5«, but, generically, non-diagonal 5Aij, this equation gives 



^ 6Xk _ SkjjSAjj _ 1 ^ / 1 l\ SAjjSAjj ^ 

k ^ - , A, A, - \A, -W Aj-Aj 

which proves (||]). 

Since matrix e is assumed to be a function of A, it can be, in fact, treated as a function 
of eigenvalues A«. Then, we use actually only 

a) the concrete form of the normalization 

b) the fact that Z[T(Aj)] is a complicated function, i.e. we should differentiate it as de- 
pending on {Aj} only through time variables. After that, ( |T8"D can be rewritten in the 



| trA 3/2 

C{VX)Z{T} 

1 



d 2 



OK 2 



d 2 Z 



Z ^ n dT a dT b 

a,b>0 a 



C(VA)ei trA3/2 Z{T} 
dT a dT b , 



<9A,- <9A,- 



56) 
^7) 



-E — 

n>0 n 



ST m^J^L^oV n ,dT n dlogC 

2>( A ^ A „ fl A.. + 2 E e ( A * 



'•J 



dAydAji 



+ 



+ 



E^)(^ T f|)frA 8 / a ) 8 -x;MA i )+ 



<9A; V3 



+ E^ 



^ ) trA 3 / 2 ) + 



/2\ <9trA 3 / 2 d/o^C 



with trA 3 / 2 = Efc A^ /2 and 



(89) 
(90) 
(91) 
(92) 



-1/2 



(93) 



The calculation of all the quantities in (pf\) - (|92|) is just an exercise of taking deriva- 
tives, using (84), all necessary details can be found in [Tj|. Careful calculation shows 
that all the terms after taking the derivatives contain only negative powers of vAi and 
can be "arbsorbed" in times. The result is: 



_2 fr A 3 / 2 

e 3 



C(VX)Z{T} [ l dA 2 



d 2 



c(VX) 



ItrA 3 / 2 



i E tr(e p A— 2 ) 

n> — 1 



fc><5n+l,0 



Z{T} 

d 2 



1 v 

a>0,f)>0 



1 



1 



m2 



^ 1 rrtrr* 



or just the set of Virasoro constraints for the case p = 2. 

In the case of generic p the analog of the derivation actually involves the same steps: 
- Represent J- [A] as 



J*»[A]=g p [A]Zto(T n ), (95) 

with 



»[A] = ^jjj nrw^«i = ^g* niAT^e*^]. (g6 ) 

- Substitute this F^[A] into (|75|) , which in the particular case of = -, 

looks like 



{Tr e(A)[( M^ )P ~ A ^9 P [MZ {p} (T n ) = 0. (97) 
d i Z 

Higher-order derivatives, { , they are defined with the help of relations like (|5J). 
- Perform the shift of variables 



T n — > T„ — T n — —S np+ % (98) 
n 

(this procedure doesn't change the derivatives). 

- After all these substitutions the l.h.s. of eq.(^) acquires the form of an infinite series 
where every item is a product of Tr[e(M)M~ k ] and a linear combination of generators of 
Wp-algebra, acting on Z^(T n ). In the case of p = 3 this equation looks like 

llY [e(M)M- 3 M- 3 "W£ ) + 
+9 £ M- 3 -V3 { E(3 * - 2)f 3 ,- 2 VV£U + E ^ W ^ 3 } + (99) 
+9E M-3n- 2 /3 | S(3A . _ 2 )T 3fc _ 2 wS_3 fc + E 



36-3 



Z& = o, 



- If N = oo all the quantities Tre(M)M~ k with given k but varying e(M) become 
independent, and (|9T| ) may be said to give H^-constraints. The exact proof for the case 



4 Integrability of GKM 



The purpose of this section is to prove that the solution to the continuum Virasoro- and 
W- constraints found above is nothing but a particular solution to the integrable KP 
system. Namely: 

(A) The partition function Z^[M] (|76|) , if considered as a function of time- variables 



T k = -Tr M~ k , k>l; (100) 

K 

is a KP r-function for any value of iV and any potential V"[X]. 

(B) As soon as V[X] is homogeneous polynomial of degree p+ 1, Z { / } [M] = zj? } [M] 
is in fact a r-function of p-reduced KP hierarchy |[20|| . [] 



In order to prove these statements, first, we rewrite ( |T6f) in terms of determinant 
formula 

ZP[M\ = det ^W „ = l,...,iV. (101) 

Then, we show that any KP r-function in the Miwa parameterization does have the same 
determinant form. | 

The main thing which distinguishes matrix models from the point of view of solutions 
to the KP-hierarchy is that the set of functions {0«(/i)} in ( |101[ ) is not arbitrary. Moreover, 
this whole infinite set of functions is expressed in terms of a single potential V[X] (i.e. 
instead of arbitrary matrix Aij in <pi(fi) = J^Aijfi^ we have here only a vector V$ or 
V[lA = 2 Vit^ 1 )- This is the origin of £_i and other W- constraints (which in the context 
of KP-hierarchy may be considered as implications of All these constraints are in 

fact contained in the Ward identity (|75|). 



7 Moreover, actually, — — — = 0. 

O-L np 



As a check of self-consistency, it can be proven (see for example ]18[) that any determinant formula 



(101) with any set of functions {<pi(n)} satisfies the Hirota difference bilinear equation. 



4.1 Integrability from the determinant formula 



We begin from evaluation of the integral: 

^ } [A] =JDX e -^V(X)-TrAX]_ (1Q2) 

The integral over the "angle" U(N) -matrices can be easily taken ^2] and if eigenvalues 
of X and A are denoted by {x{\ and {\{\ respectively, this integral can be rewritten as 



A(A) 



V N 

1] / d Xl e- v ^ )+x ^ 

.2 = 1 J 



A(X). 



A(X) and A(A) are Van-der-Monde determinants, e.g. A.(X) — Yli>j( x i ~~ x j) 
The r.h.s. of ( |103| ) can be rewritten as 





with 



Note that 



If we recall that 



A- 1 (A)A(^-)n J dx ie - y ^ )+x ^ 



A- 1 (A)det fo) F i (A j/ 



F i+1 {\) = Jdx x l e~ v ^ +Xx = (^)^i(A) 



Fx(A) = 4S[A] 



(103) 



(104) 



(105) 



(106) 



A = V'(M) 

and denote the eigenvalues of M through {//j} , then: 



(107) 



with 



(108) 



Proceed now to the normalization flTF). Indeed, it is given by the Gaussian integral: 



C^[V\M]- X = J DX e~ u ^ M ' x \ (110) 

Making use of ?7(iV)-invariance of Haar measure dX one can easily diagonalize M. Of 
course, this does not imply any integration over angular variables and provide no factors 
like Apf). Then for evaluation of fllToD it remains to use the obvious rule of Gaussian 
integration, 

r ^ N N 

/ DX e ~^i u ^ x " ~ l[U^ 1/2 (111) 

J id 

and substitute the explicit expression for C/j,-(M). If potential is represented as a formal 



V(X)=Y— X n , (112) 

n 



we have 



and 



U 2 (M, X) = J2 v n+i \ E TrM a XM b X 

n=0 I a+b=n— 1 



n=0 la+6=n-l I n=0 W W ^ Mi 



Coming back to (|76|) , we conclude that 



4 y} [M] = e Tr [ y ( M )- w '( M ^W[V|M]^[V'(M)] ~ 



sfri) = [y"(M)] 1 / 2 e F ^-^ y '^ (114) 
The product of s-factors at the r./t.s. of ( |113| ) can be absorbed into ^-functions: 



AV\ r ,^ det$i(u 7 ) 



where 



(116) 



where the asymptotic is crucial for the determinant (|1 1 5|) to be a solution to the KP 



hierarchy in the sense of pO 



The Kac-Schwarz operator [p3j 24j| . From eqs. (|109|) , (|114|) and (|H6|) one can deduce 
that $j(yu) can be derived from the basic function by the relation 



$.(^) = [V"(fi)} 1/2 J x l - x e- y( - x) + xV '^dx = Apyi^in) 
where l{y}(/z) is the first-order differential operator 



(117) 



i 



d 



E^(aOP 9^ [^"(/i)] 1/2 



In the particular case of V(x) 



V"(fi) dfx 



rP +i 



p+1 



Id p — 1 

+ /i 



118) 



p/jP" 1 dfi ' r 2pfi p 

coincides (up to the scale transformation of fi and A[ p }(/x)) with the operator which 



determines the finite dimensional subspace of the Grassmannian in ref. [f2I| We emphasize 
that the property 



$ m (/i) = V}(aO^) {Fi+iW = -§x m)) 



(119) 



is exactly the thing which distinguishes partition functions of GKM from the expression 
for generic r-function in Miwa's coordinates, 



Tn m ~ A(M) 



(120) 



with arbitrary sets of functions <f>i(fi). In the next section we demonstrate that the quan- 
tity (11201) is exactly a KP r-function in Miwa coordinates. 



4.2 KP r-function in Miwa parameterization 

A generic KP r-function is a correlator of a special form [|15j : 



r G {T n } = (0| : e£ T " J " : G\0) (121) 

with 

J(z) = i)(z)ip(z); G = : exp Q mn ^m^n ■ (122) 

in the theory of free 2-dimensional fermionic fields ip(z), ip(z) with the action / ipdip. The 
vacuum states are defined by conditions 

ip n \0) = n < , $ n |0) = n > (123) 

where ip{z) = Ez i^" ^ 1/2 , i>{?) = Ez ^nZ~ n-1 dz 1 / 2 . 

The crucial restriction on the form of the correlator, implied by ( |122j ) is that the 
operator : TnJ " : G is Gaussian exponential, so that the insertion of this operator 
may be considered just as a modification of (ipip) propagator, and the Wick theorem is 
applicable. Namely, the correlators 

(0|nfe)V>(Ai)G|0) (124) 

i 

for any relevant G are expressed through the pair correlators of the same form: 

(124) = det fo - ) (0|^(/i i )^(A,)G'|0) (125) 

The simplest way to understand what happens to the operator e^ TnJn after the sub- 
stitution of ( |10 (JQ is to use the free- boson representation of the current J(z) = dcp(z). 
Then T,T n J n = X) — ~ir^W = XM^)> and 

:e EM^i).- J rjjie^:. (126) 

In fermionic representation it is better to start from 



instead of ( |100| ). Then 



= eET " J ' 1 := n Tl ^r/rr _ r , n^M^) ■ (128) 

In order to come back to ( |100|) it is necessary to shift all /Vs to infinity. This may be 
expressed by saying that the left vacuum is substituted by 

(N\ ~ (Ol^oo^oo)../!/^- 1 ^). 
The r-function now can be represented in the form: 

t%[M] = (0| : e £ T " J " : G\0) = A(M)" 1 (AT| J] : : G\0) = 

n-(A-M) ~* (129) 

« rii>j c/x< - ^ ) rii>j c^t - » 

applying the Wick's theorem ( |124j) , (|125|) and taking the limit fa — > oo we obtain: 

T =[M] = (13 0) 

with functions 



<M/i) ~ (OI^M^CaiJGIO) -> ^(1 + 0{-)). (131) 



Thus, we proved that KP r-function in Miwa coordinates ( |100| ) has exactly the deter- 
minant form (ITH), or, put differently, ( |101| ) is a r-function of KP hierarchy. Below we 



will discuss how from generic point of Grassmannian described by G = exp Aijipiipj or 
infinite matrix with two indices (oo 2 ) one can restrict it to a solution, determined only by 
one function (oo) or two functions (2 x oo).. 



5 Universal £_i-constraint and string equation 

Let us return to the question of specifying particular "stringy" solutions to the KP hier- 
archy. It turns out that using integrability it is enough to prove only the so-called string 
equation or £i-constraint, all other recursion relation follow from these two statements 



It is well known that £_i-constraint is closely related to the action of operator 



d Id 

Tr^— = Tr — -— . (132) 

dA tr V"{M)dM tr v ; 



Therefore it is natural to examine, how this operator acts on 



siji) = (y^jV-^'W, (134) 



$,(/i) = Fi(X) = (d/dXy-'F^X), X = V'(fi). 
First of all, if is considered as a function of T-variables, 

1 rp d {v} 1 dlogZW 

ZW dA tr " V"{M)M n + l1 dT n ' 1 j 

On the other hand, if we apply (|132|) to explicit formula ( |133|) , we obtain: 



1 r 2V-?_zW 



Zi y } dA 



tr 



.Tr M + i £ , ' , Z^'-^) + T r #- log det F,( V 



^XjV'i^V"^) in -fij OA 



tr 



(136) 



Below, we will argue that 

_L £ m z m = _ JL log Z {v} + TrM _ Tr J_ log det F . (A . } . (137) 

can be used in order to suggest the formula for the universal operator C^Q . 
Here 



i 



d 



1 \- 



„ , L V"(M)M^ 1 dT n 
V'iin) - V'\to) d 



+ 



2^V"(to)V"(to) 



^ - to 



(138) 



and this turns into more common expression when V(X) = X p+1 /(p + 1) (note that the 
items with i = j are included into the sum at the r.h.s. in (|138|) ). 



So, in order to prove the £^-constraint, one should prove that the r.h.s. of (|137|) 
vanishes, i.e. 

log ZSP = TrM - Tr— - log det Fi(Xj), (139) 



m ° " dA tr 

This is possible to prove only if we remember that = . In this case the l.h.s. 

may be represented as residue of the ratio 

T^ } {T n + fi- n /n) d {V} (rr ^ / 1/m x 

res^ ^y-— = log r x N f (T n ). (140) 

T N ( T n) 01 1 

However, if expressed through Miwa coordinates, the r-function in the numerator is given 
by the same formula with one extra parameter fi , i.e. is in fact equal to rj^ 1 . This idea 
is almost enough to deduce ( |139[ ). Let us begin with an illustrative example of N = 1. 
Then (A = V'Qj,)) 



r[ v \T n ) = rPi^) = e^^ v '^[V"(» 1 ))V 2 F(\ 1 ), 



Tl v \T n + ^ n /n)=rP[^^} = 

= ^y^v'M^)-^) \y"MV"W ,% [ F ( Al )aF(A) /dX - F{X)dF(X 1 )/dX 1 

ji — Hi 

= —rPW] ■ [-dlogF(X 1 )/dX 1 + dlogF(X)/dX]. 

H — Hi 

The function 



(141) 



r V"" 
F(X) = J dx e - y{x)+Xx ~ e v M-v v 'M[v , '(jj l )]- 1 / 3 {l + 0(_1__)}. (142) 

If V(fi) grows as fi n when fi oo , then V""/(V") 2 ~ fi~ n , and for our purposes it is 
enough to have n = p+ l>l,so that in the braces at the r.h.s. stands {1 + o(l//i)}(/i • 
o(/j,) — > as — ► oo). Then numerator at the r.h.s. of Ql41| ) is ~ 1 + o(l//x), while the 
second item in square brackets behaves as dlogF(X)/dX ~ /i(l + o(l///)). Combining all 



9 \ogrP = res, { 1 + [-dlogF(X 1 )/d\ 1 + /x(l + o(l//x))] 



= fii-dlogF(X 1 )/d\ 1 . (143) 

i.e. ( |139| ) is proved for the particular case of N — 1. 

The proof is literally the same for any iV , we will omit here the details which can be 
found in ||18|| . After a simple but cumbersome calculation one gets 



9 i {V} 

- log r x N 



1 + 0(1/0) 



™>{ u^£f {[1+0(1/ ^ ] - 

^[Tr^|- log det Fj(Aj)] • [1 + 0(1/ //)] 

iV a 

^ /i ^ Tr ^ 1 ° g deti ^( A i)" ( 144 ) 



which completes the proof of eq.(3.13) and thus of the universal C_x -constraint. 

In the particular case of monomial potential V = V p = XP ^ ( |138| ) turns into more 



p+i 

common form M. |5f : 



+i £ aT.&T,, - A (145) 

4P a+b=p U - L l 



6 GKM versus Toda theory and discrete models 

Now, first, without any special reference to GKM there exists an explicit relation between 
KP-like (in Miwa variables), 

TKp[Tt] = (M6) 



T~N[T-k,Tk] — det Hi + N,j+N[T-k,Tk], 

1-3 

representations of r-functions, where 



(147) 



A( f i) = H^ t -fi j ), (148) 

i>j 



^)=^-\l + 0{-)l (149) 



T * = ^E/^fc>0, (150) 



3H l3 /dT k = H iti _ k , j>k>0, (151) 

and 

dHi./dT^ = Hi- kd , i>k>0. (152) 
Relation between ( |146| ) and ( 147 ) is formulated in terms of the Schur polynomials, 



which are defined by: 

P[z\T k ] = exp{£ T k z k } = £ z k P k [T], (153) 

k>0 

e.g. P_ n = for any n > 0; P [T] = 1; P^T] = T\; P 2 [T] = T 2 + §7?; P 3 [T] = 
T 3 + T 2 Ti + |T X 3 etc. The crucial property of Schur polynomials is: 

dP k /dT n = P k . n (154) 

(this is just because dP/dT k = z k P). This feature allows one to express all the depen- 
dence on time- variables of Hij[T], which satisfies eqs.( |151| ) and ( |152|) , through the Schur 
polynomials: 

ifypZLp, T p ) = P- k [T- p }H kl P l+J [T p } } (155) 

k<i 
l>-3 



Let us begin our consideration from the case, when all N = T_ k = 0, then look what 
happens if iV > f\, and introduce TL^-variables. 

Given the system of basic vectors for i > 0, we put by definition 



0,T k 



z^O 



[z)z- j V[z\T k }dz, i > 0. 



(156) 



The integration contour is around zero and it can be deformed to encircle infinity and the 
singularities of V[z], if any If we just substitute the definition ( |153| ) of V[z] into ( |156| ). 
we get (fjM) with P k _i[T- m = 0] = 5 ki and 



kl 



z^O 



4> k (z)z l dz. 



(157) 



In order to prove the identity between (|146|) and ( tL47|) let us substitute the Miwa trans- 
formation (jsg) to (p3|) 



V[z\T k ] 



det M 



n 



det(M-Iz) V (fa - z) 
- fj,j). Now t 
from the poles of V[z\T k ] at the points 



i J k 



E 



(~) fc A fc (/i) _ 
(z - /i fc ) A(/i) ' 



where A fc (/i) = II (/ij — //,). Now the integral ( |156| ) picks up contributions only 



HiAT- 



0,71 



z^O 



[z)z->r[z\T h ]dz = ^U-nM^ 



The sum at the r.h.s. has a form of a matrix product and we conclude that 



det Hy = det • J ] 



A(//)' 



: A fc (/i) 



• det 



/4 



The last determinant at the r.h.s. is equal to A(l//x) ~ A(/x) • fife 



. Note also that 



n 



k L 



A fc (/i) 

A(//) 



A(/i) 2 , and taking all this together, we see that there exists the equality 



det 0^) 
d6t ^ = AM ' 



158) 



as required. 



3 discussion of N < can be found in 



Proceed now to introducing of zero- and negative-time variables. The zero-time n 
arises just as the simultaneous shifts of indices i and j of ify : — > H i+N j +N . We can 
write: 

H i+N>j+N [0,T k } = <f 4 N} (z)z- j V[z]dz (159) 

with 

4 N \z) = z~ N ^{z) (160) 



This exhausts the problem of restoring the iV-dependence for positive integer values of 
N. 

As for negative-times, as soon as H k i is defined, they are introduced with the help of 
and 



Hi + Nj+N[T-k,Tk] = Pj-k[T-p]H k+N ^ +N \0, T t ] 



k<i 



with 



= / <pf T - k ' N} (z)z^V[-\T^ k }V[z\T k ]dz, 

Jz^O Z 

= z~ N ex P \ -Y,T_ k z- k \ Y,Pk[T-i]<t> i+ N-k{z) 



(161) 



(162) 



fc>0 

The role of the exponential prefactor in ( |162j ) is to guarantee the proper asymptotic 
behaviour 

4 T - k ' N \z) = z i - l {l + 0{-)}. (163) 

z 

The important reduction from Toda lattice is Toda chain (see, for example, p0|). It can 
be easily written both in terms of element G in the fermionic language ([G, J k + J- k ] = 
- see also the comments above, and in determinant form. Latter one merely implies the 
symmetry property: 



where A is shift matrix A^- = This condition leads to r-function of Toda chain 

hierarchy (proper rescaled by exponential of bilinear form of times) which depends only 
on the sum of positive and negative times t k = \{T k + T_ k ), but not on their difference 
(one can consider this as defining property of Toda chain hierarchy). Let us remark that 



one possible solution to constraint (|164|) is matrix = Tii-j. We can combine both 
reductions to reproduce forced Toda chain hierarchy. In this case one can easily transform 
Hi~j to matrix Tii+j by permutations of columns what does not effect to the determinant. 
This matrix just corresponds to one- matrix model case [12], Thus, we consider again 



the determinant of size N x N, which now can be represent in the form: 

t n = det d i+i H (165) 

where d = d/dt±, dH/dt k = d k H. Like the Toda lattice case, the forced Toda chain is 
unambiguously continuable to negative values of zero-time. 

Now, we can easily introduce zero- and negative-time variables into GKM in such a 
way that its partition function becomes a r-function of the Toda lattice hierarchy. The 
relevant set of functions {4>i{fi)} - the point in Grassmannian for GKM - is given by the 
following integral formula: 



(166) 



EE 8(JJL) J dx x ^ e -VW+xVto = 

Dependence of N and T_ k is now introduced by the rule P|: 











X 


" J 






exp < 









l>0 




yJV"((x)e v M » v M f ^ ^ jV+i _ 1 ^_ F(a:)+ „ y , (M) _ _ ^| _ (167) 



f dx x N+i-i e -v( x )+ x v'M exp J ^T. 



10 Let us point out that the exponential of negative powers in normalization does not essentially ef- 
fect to the KP r-function as it reduces to trivial exponential of bilinear form of times in front of r- 
function and corresponds to the freedom in its definition. Indeed, t ~ det {exp[^_! fe a k zj ]<f>i{zj)} ~ 



where 



V(X) = V(X) -N\ogX-J2 t- k X~ k (168) 

fc>0 

with 

The original potential V it can be identified with V + . From (|167|) we immediately 
conclude that the partition function of GKM, involving zero- and negative-times (and 
automatically being a Toda lattice r-function), is just 

f DY „-TrV(X)+TrV'(M)X 

Z m \M\ = e TrV(M)-TrMV;(M)U^ (16g) 

{^> L J J dX e -TrU+, 2 (X,M) V ; 

Since we devote this section to discussion of Toda lattice hierarchies in the context of 
matrix models, we can not avoid touching the main conclusion of that all the discrete 



matrix models do correspond to particular cases of Toda hierarchies. In the simplest 
case of Hermitean one-matrix model one gets a Toda chain, other multi-matrix models 
correspond to other reductions of the Toda lattice hierarchy. Moreover, all discrete matrix 
models fall into the class of forced hierarchies JT3| . 



The idea is to perform a Miwa transformation of times with fixed, so that Hi/s 
become averages of polynomial functions of X. Then det Hij may be transformed with 
the help of orthogonal polynomials technique. 

The main result of all these calculations is that partition functions arise in the form 



( 146 ), with <fii(z) being proportional to orthogonal polynomials. 
Indeed, using 

1 1 R 

T k = -TrA- k = - £ K k + tk (170) 

(k > 0) □ (note that N - the size of matrix A has nothing to do with N - the size of 
matrices M, being integrated over in ([18]) ) the partition functions of discrete models 
{tn} acquire the form of KP r-function ( |146| ). E.g., for the Hermitean one-matrix model 



one has 



11 and/or 

1 1 A 

-TrM- k = ■■- k 
k k 

i=l 



T-i, = T, : r Tr\J ' = T 



r N {t) = (Nl)- 1 ! Y[dm,A 2 (m)ex P {-J2tM} 

i i,k 



(N\)- 1 I n^ i A 2 (m)e- v ^n( 1 ~ T~) 

J _• _ A„ 



(iV!)- l nA- iV /n^e-^A( 



A(m, A) 
'^(AT 



(A^!)~ 1 nA^ A 'A- 1 (A) /n^e-^ (mi) x 



x det Pi-Arrij) det 

JVx^V (N+N)x(N+N) 



P&K) : P 



N + b-l{ m j, 



(171) 



where z, j = 1, . . . , iV; a, 6 = 1, . . . , N; and {Pj(m)} are corresponding orthogonal polyno- 
mials with respect to deformed measure e~ v dm : 



< P^Pj >= J P i (m)P j (m)e-^ (m) rfm = <5 -e^ w , 



(172) 



so that 



P(m) = m ! + 0(m^ 1 ). 
Computing determinants in (171) and using orthogonality condition (|172j ) one obtains 



T N [\\i\ = H\- N A-\\) det PAr + a-i(A b )II e 



<Pi(f) 



NxN 



n 



det (o6) 0W(A 6 ) det (ab) 0W(A 6 ) 
- = Tat [t\ x - 



(173) 



A(A) JVLJ A (A) 

i.e. the r-function of the discrete Hermitean one-matrix model acquires the form of 

eq. fljgp with 



<P[ N \\) = \- N P N+a - 1 (\) (174) 

(|173 ) is natural representation for all discrete matrix models. 

Now, in order to be representable as GKM the above formulas still need to arise in 
a somewhat specific form. Namely, components of the vector {0«(/i)} should possess a 



Since in the study of discrete matrix models <pi(fi) arise as the orthogonal polynomials 
P n+i (fi), what is necessary is a kind of integral representation of these polynomials, with 
z-dependence coming only from the x J ~ 1 -factor in the integrand. It is an interesting 
problem to find out such kind of representation for various discrete models, but it is 
easily available only whenever orthogonal polynomials are associated with the Gaussian 
measure: the relevant Hermit polynomials are known to possess integral representation, 
which is exactly of the form which we need. 

The final statement we have got here is that Hermitean one-matrix model with the 
matrix size N is equivalent to GKM with V(X) = X 2 /2 - NlogX. 

Indeed, we have proven the remarkable explicit identity between Gaussian matrix 
integrals, 



JDM NxN det(I-M/A)e- TrM2 / 2 _ J M M det(Z - iX / X) N e~ TrX2 ' 2 

JDM NxN e-^MV2 !DXMe -TrX V 2 y ) 

Note that the size of matrix in the l.h.s. is N x N and in the r.h.s. is N x N, and these 
parameters are absolutely independent. This identity is indeed true for any N and N, 
as follows from the proof given above and integral representation of Hermit polynomials 
4>j{(i) = Hej{ifi). 



7 Double-scaling limit 

Now, we are going directly to the discussion of the connection between discrete and 
continuum theories. Indeed, it turns out [^7]] that this connection at the language of free 
scalar fields is nothing but a change of spectral parameter 



u 2 = 1 + az (176) 
Then the continuous Virasoro constraints (|5]) which are modes of the stress tensor 



can be deduced [27J from analogous constraints in Hermitian one-matrix model by taking 



the continuum limit. The procedure is as follows. 

The partition function of Hermitian one-matrix model (|2|) satisfies the discrete Virasoro 
constraints (|5[). In order to obtain the continuum constraints Q78|) one has to consider a 
reduction of model fll8D to the pure even potential t 2 k+i = 0. 

Let us denote by the r^ d the partition function of the reduced matrix model 



T v N ed {t 2 k} = J DMexpTr^^M 2fc (178) 



k=0 

and consider the following change of the time variables 



9m- £ Lf \i\ (179) 

n >m [n-m)\r m+o 



2, 

(g m = mtim and this expression can be used also for the zero discrete time g = N that 
plays the role of the dimension of matrices in the one-matrix model), deduced from the 
following prescription. Take the free scalar field with periodic boundary conditions for 
p = 2) 

dcp(u) = £ g k u^ + £ ^-u- 2k -\ (180) 

fc>0 k>l ° t2k 

and analogous scalar field with antiperiodic boundary conditions: 

&>{z) =Y,(( k + l) W**"* + ^—z- k -i) • (181) 
fc>o\ v 2J dT 2k+1 j 



Then the equation 

~d$(z)f = a— 

generates the correct transformation rules ( |179| ), ( |184|) and gives rise to the expression for 



■^d${z)f = az ^dv(u)T md , u 2 = l + az (182) 



A nm ( P-87j ). Taking the square of the both sides of the identity ( |182j ) 

4 _ 1 



^T{z)f = z^T( M )r red , (183) 



T T L 



one can obtain after simple calculations that the relation ( |188|) is valid. 
Derivatives with respect to t<ik transform as 



dUh (h - -n - 1 MP (n M ^ (9To 



where the auxiliary continuum times T2 n +i are connected with "true" Kazakov continuum 
times T 2n+1 via 



k 



185) 



•jfe + l/2- 2( *- 1)+1 ' 
and coincide with T2 n +i in the double-scaling limit when a — > 0. 

Let us rescale the partition function of the reduced one-matrix model by exponent of 
quadratic form of the auxiliary times T 2n+ i 



1 



f — exp ( — — ^2 A mn T 2m+ iT 2n+ i ] t n 

m,n>0 



.red 



with 



A, 



r(Wf)r(m + 



hmg-n-m-l 



2r 2 (j^J n\m\(n + m + l)(n + m + 2) 

Then a direct though tedious calculation |27| demonstrates that the relation 



(186) 



(187) 



n+1 



^red_red 



ECf-i)"" 1 ^ 

Z— / n +J- V / ^-red 



p=0 



is valid, where 



- red 
J 2n 



k=0 







2 k 



dt 



+ E 



d 2 



(188) 



2(fc+n) < fc < n dt 2 kdt 2 (n-k) 



(189) 



and 



£_i= Wfc + ^T 2fc+1 
fc>i v z/ 



fc>0 v 



5 



+ 



ar 2(fc _ 1)+1 " i6G' 
9 



4= E ^ + ^) T2fe+1 ^i^ 



fc>0 

+ E 



5 



5 



Here Cl 



0<fc<n-l (n-fc-l)+l 

are binomial coefficients. 



16a r 



n > 1. 



(190) 



J ri p\(n-p)\ 

These Virasoro generators differ from ([f^) by terms which are singular in the limit 
a — > 0. At the same time U 2 d r red at the r.h.s. of (|188|) do not need to vanish, since 

= L 2p r t2k+l=0 = U 2 fr^ + £ — jf r . (191) 



It was shown in P7[] that these two origins of difference between (|78| ) and ( |190| ) actually 
cancel each other, provided eq.( |188| ) is rewritten in terms of the square root y/f rather 
than f itself: 



n+1 

a n ( 
p=0 



m+l-p -^2p r 



;i + 0(a)). 



(192) 



The proof of this cancelation, as given in ||27|| , is not too much simple and makes use of 
integrable equations for r. 

Now we will use the demonstrated above fact, that the discrete Hermitean one- matrix 
model is equivalent to GKM with V(X) = X 2 /2—NlogX, and also that its double-scaling 
continuum limit is described by GKM with V(X) = X 3 /3. Thus, we should conclude that 



lim Z ( y, 

d.s. N-+00 *- v J 



<{vy 



(193) 



This relation should certainly be understandable just in terms of GKM itself. 



Let us recall that double-scaling continuum limit for the model of interest implies that 

only even times /^fc — — Tr^-r should remain non-zero, while all odd times t 2 fc+i — 0- 

2k A 2k 

This obviously implies that the matrix M should be of the block form: 



A 



(194) 



' M 
v —M 

and, therefore, the matrix integration variable is also naturally decomposed into block 
form: 



X 




(195) 



Then 



Z{V=X 2 /2-N\ogX} ~ 

1 (196) 

= J DXDyD 2 z det(xy- z^zyf e- T ' r ^ z \ 2+x2 ' 2+y2 ' 2 - MX+My \ 

To take the limit — > oo, one should assume certain scaling behaviour of X, y and Z. 
Moreover, the notion of double-scaling limit implies a specific fine tuning of this scaling 



X = a(i/3I + x), 

y = a(-ipi + y), (197) 
Z = a(, 

with some large real a and /3. 

As for behaviour of the matrix of Miwa's parameters, it should be dictated by the 
change of the spectral parameter (|176|) , which dictates the right Kazakov change of vari- 
ables. So, we will take the anzatz 

M 2 = A + Bm (198) 
which under particular scaling behaviour A/B — > oo could be linearized giving rise to 

M = a~ l (i^I + m) (199) 

with 

7 



A l/2 = 1 

a 

B 1 



a 

If expressed through these variables, the action becomes: 



Tr{\z\ 2 + x 2 /2 + y 2 /2 - mx + My - N\og(xy - z^zy)} 



(200) 



—Tr{{if3I + x) 2 + —Tr{if3I - y) 2 + 7 2 |C| 2 } - Tr(ial + m)(2i(3I + x - y)- 



-NTr log /?V{1 - ^ + P - + o(l//?))} = 

(201) 

= [2a/? - /?V - 2iV log (3j]Tr I - 2i/3 Tr m+ (A) 
+*(/3 7 2 -a + ^)(Trx-Trj/) + ^( 7 2 - —){Tr x 2 + Tr y 2 )+ (B) 



Hl 2 ~pTr\t\ 2 - (C) 
iN 

—Trmx + Tr my + ^^Tr(x 3 — y 3 )+ (-D) 



We want to adjust the scaling behaviour of a, (3 and 7 in such a way that only the terms 
in the line (D) survive. This goal is achieved in several steps. 

The line (A) describes normalization of functional integral, it does not contain x and 
y. Thus, it is not of interest for us at the moment. 

Two terms in the line (B) are eliminated by adjustment of a and 7: 

9 N 2N , , 

t 2 = ^.« = T - <*»> 

As we shall see soon, j 2 = N//3 2 is large in the limit of n — > 00 . Thus, the term (C) 
implies that the fluctuations of ^-field are severely suppressed, and this is what makes 
the terms of the second type in the line (E) negligible. More general, this is the reason 
for the integral Zcy-, to split into a product of two independent integrals leading to the 
square of partition function in the limit n — > 00 (this splitting is evident as, if Z can be 

I x z\ 

neglected, the only mixing term logdet turns into \ogXy = log A? + log y). 

\z y) 

Thus, we remain with a single free parameter (3 which can be adjusted so that 

3 3 

j- -> const as iV -»• 00 (i.e. (3 ~ N 1/3 ), (203) 
making the terms in the last line (E) vanishing and the third term in the line (D) finite. 



Let us finish the discussion of double-scaling limit, making some remarks on ref. [pq] . 
It is claimed there that in order to get the Kontsevich model from a discrete Hermitean 
one, one should not necessarily care of the reduction to even times ( |194| ) and the particular 



Kazakov change of variables, inspired by (|176|) , (|198|) ; instead it is just enough to take 



(|199 ) and make the rescalings similiar to what we done above. However, in such case it is 



not clear what should be instead of Kazakov change of variables, and what is indeed the 
integrable hierarchy we get in continuum limit. As for the first question, one might hope 
that the related change of variables in in the class of "allowed redefinitions" of Kazakov 
times, similiar to those we used in ( p.85| ), though the second question, certainly deserves 
further investigation. 



8 GKM as a solution to topological (p, 1) models 



For various choices of the potential V{X) the model (|7q) formally reproduces various 
(p, g)-series: the potential V(X) = can be associated with the entire set of (p, q)- 
minimal string models with all possible g's. In order to specify q one needs to make a 
special choice of T- variables: all = 0, except for T\ and T p+q (the symmetry between 
p and q is implicit in this formulation). 

Let us briefly discuss two simple examples. First, we will fix p = 2, i.e. the case of the 
Kd V reduction to the KP hierarchy. The second number q should be coprime to p, thus 
we have here q = 2m — 1. 

Now the string equation is of the form 

' C-xTKdv = \ E kT klJ— ^r K dv + % = (204) 



fc odd 

or taking the ^-derivative, one gets 



E kT k — log r KdV + T 1 = (205) 

fc >i Olk-2<Jll 

k odd 



or using the formula for the Gelfand-Dikii polynomials 

d 1 



we have 



\ogr KdV = |L J ^ = Tl m [u\ (206) 



E (2m + l)T 2m+1 K m [u\ = (207) 



m>0 



Now, we should use the "axiomatics" of Q how to extract the concrete (2,2m — 1) 



solutions from ( |205|) , (|207| ). The very simple example is m — 1: where 



3T 3 ^logr KdV + T 1 = (20£ 



using that 

& 



dT{ 

the solution to the KdV-ennation is 



or fixing T 3 , one gets 



(209) 



F = log r~Tf (210) 
which is a well known fact from the c = —2 theory coupled to gravity where 

(P 3 ) = 1 

with P being the puncture operator P = cce^. This is the case of topological gravity. 
Less simple example is with m = 2, where Q207| ) and the fact that 

T> 2 i // 

gives the Painleve equation 

u 2 + u" = Ti (211) 

This is the case of pure gravity where the solution is actually much more complicated 
than in the previous case. 

From this point of view, the presense of all (p, q) solutions in GKM is a rather formal 
consideration. For the potential V(X) = the partition function Z[V\Tk\ = TyfTfc] = 
r p [T k ] satisfies the string equation which looks like 

P-1 oo v d 

^k(p-k)T k T p „ k + ^(p + k)(T p+k -^6 kjl ) — log r p [T}=0 (212) 
k=i k=i P + 1 01 k 

i.e. r-function is defined with all Miwa times ( |100| ) around zero values (in 1/M decom- 
position like in original Kontsevich model) with the only exception - T p+1 is shifted what 
corresponds obviously to (p, 1) model. Thus, we see that the matrix integral gives an 
explicit solution to (p, 1) string models which must be nothing but particular topological 
matter coupled to topological gravity. 

Of course, we still have an opportunity for analytic continuation in string equation, 
using the definition of Miwa's times (|100|) . We have to satisfy the following conditions: 



T 2 = 



T P 



p+1 

Tpj^q tp-\-q fixed 

T p + q +l = 

(213) 

which is a system of equations on the Miwa parameters {/ij}, % — 1, N. So, to do this 
analytic continuation one has to decompose the whole set 

{&} = {£«} © K} 
1 1 N 1 1 N ' 

T k = -TrM~ k = lj2^J k = T.J2 C k + tE ^7* = + T' k (214) 

3=1 3=1 

into "classical" and "quantum" parts respectively. In principle it is clear that we have 
now to solve the equations 



and this can be done adjusting a certain block form of the matrix M fig , |30| . However, 
in such a way we can only vanish several first times, the rest ones can be vanished only 
adjusting correct behaviour in the limit iV — ► oo. The most elegant way to do this is 
to use the formula 

ex P (- £ A fc T fe (c/) ) = lim (1 - ±- £ A fc lf >)* = H(l - f ) (216) 

k=l A fc=l a ?o 



and then the solution to (|215 ) will be given by K sets of roots of the equation 



£ A fc T fc (c ° — K = t p+q X p+q ^-A p+1 - K = (217) 

fc=i p+1 

Obviously, the eigenvalues £ a will now depend on the size of the matrix N = (p + q)K + N' 
through explicit independence (£ a ~ ffvteH-?)) aric l we i ose one Q f the main features of 
(p, 1) theories mentioned above - trivial dependence of the size of the matrix. Now we 



can consider only matrices of infinite size and deal only with the infinite determinant 
formulas. 

That is why this way to get higher critical points is a formal one. Below we will 



present an alternative way of thinking [29, pry, p2|, connected with so-called p-times. 



Indeed, there exists a priori another integrable structure in the model ([7f]), connected 
with time variables, related to the non-trivial coefficients of the potential V. As a results, 
the cases of monomial potential V P (X) = ^^j- and arbitrary polynomial of the same 
degree (p + 1) are closely connected with each other. 

In order to demonstrate this, first, we return to the derivatives of Zqkm with respect 
to the time- variables T^. Such derivatives define nonperturbative correlators in string 
models and are of their own interest for the theory of GKM. The derivatives with respect 
to Tfc with k > p + 1 (responsible for the correlators of irrelevant operators) are not very 
easy to evaluate, things are simpler for T k with 1 < k < p, where using the obvious 
notation of average so that Zqkm = (1), we have 

dZ GKM 



= (TrM k - TrX k ) , 1 < k < p (218) 

v 



It is implied that the derivative in the l.h.s. is taken preserving the form of the potential 

V = E p+1 f x k . 



The r.h.s. of (|218| ) can be also represented as 



dZ GK M 



l Tr d -m-Tr^l\^<k 



v 



<p (219) 



which looks similar but actually is different from —^-Zqkm-i as h would be if ([219D 



does contain some corrections. The problem is that -£^Zgkm gets contributions not 
only from differentiating V(X) — V(M) in exponentials in (|76| ) but also from the term 
V'(M)(X-M) = W{M)(X-M) as well as from the pre-exponential fl77|). The corrections 
are of the two different types 

d 

— W) + "quantum corrections" (220) 
dv k 

First type should dissappear if one introduce a new "spectral parameter" 



Thus, we are led to special time variables induced by a special transformation of the 
spectral parameter \i 



f k = \TrM- k (222) 

The formula ( |220| ) also demonstrates that the {v k } are not the true time variables for a 
given arbitrary potential. Indeed, it appears that the right variables are the parameters 
{t k } being certain linear combinations of the coefficients {v k } of the potential |35|, |36[ 



k(p — k) 

Using (|223|) one can get two important formulas 



?—Res W x - k/p {ij)dn (223) 



i p+i 

/i = -5>t fc/ i fc -P, (224) 

and 

p+i 

V( f i)-fiV'(fi) = -Y / tk{i k . (225) 

— oo 

The second one implies the natural interpretation of the exponential pre-factor in eq. (fT6p 
as the standard essential singularity factor in the Baker- Akhiezer function of p-time vari- 
ables. 

Now, the direct calculation shows that 

Z[V\T k \ = r v [T k ] = 
= ew(-\Y,Mt)(T l -U)(f J -t j )y p if k -t k ] , (226) 

where 

A i; j = Res^W i/p dWi /p , (227) 

and /(yu)+ denotes the positive part of the Laurent series f(fi) = S It is also easy 
to demonstrate, that 

r p [T}^r Vp [T] (228) 

- is the r-function of p-reduction. 

Formula ( |22(j| ) means that "shifted" by flows along p-times ( f^23| ) r-function is easily 



time- variables and The change of the spectral parameter in ( [21 3| ) M — > M 



f(M) = W 1 / p (M) (and corresponding transformation of times Tjt — > 2\) is a natural step 
from the point of view of equivalent hierarchies. 

Indeed, the relation between r-functions of the equivalent hierarchies can be easily 
derived from an identical transformation: 

r(T) = ^Y[[f(^f(f) (229) 
where f(T) as function of times T has the determinant form (|101j) with the basic vectors 

m = if'ww^w)) (230) 

By a direct calculation one can show that pre-factor in eq.( [229| ) may be represented in 
the form 

^n[f(^)] 1/2 = exp f-^EA^j (231) 

where 

A l3 = Res f\X)d x fi(\). (232) 



From ( 229 ) we see that 



r(T(T)) = f(T)exp ( -Wa^ ] (233) 



Let us introduce the r-function f(T) of the p-reduced KP hierarchy defined by: 

f(T)^^|exp [$> fl) 



j 



T 



(t)= e -|E^*J (234) 



for which instead of 



we have 



_ det ^(/i,) 

rv[T] - A(//) (235) 



r p [T -t] _ det 4>i(fij) 



r p [t) A(/2) 

with the corresponding points of the Grassmannian determined by the basic vectors 



(236) 



r 



and 

iiip) = [pjl p - 1 } 1/2 exp ( - V tjjlj \ I x^e-^+^dx (238) 




respectively. Then it is easy to show that f p (T) satisfies the L_i- constraint with shifted 
KP-times in the following way 

p-1 OO Q 

HP - k)(f k - t k ){f p „ k - tp_ fe ) + J2(P + k)(f p+k - t p+k )-^logr p [f - t] = (239) 
fc=i fc=i ol k 



where ti defined by ( |223| ) are identically equal to zero for i > p + 2. 

The formulas fl226|J239p demonstrate at least two things. First, the partition function 
in the case of deformed monomial potential (= polynomial of the same degree) is ex- 
pressed through the equivalent solution (in the sense [|39], (40]) of the same p-reduced KP 
hierarchy, second - not only t p+ i but all t k with k < p + 1 are not equal to zero in the 
deformed situation. We will call such theories as topologically deformed (p, 1) models (in 
contrast to pure (p, 1) models given by monomial potentials V P (X)), the deformation is 
"topological" in the sense that it preserves all the features of topological models. More- 
over, this "topological" deformation preserves almost all features of 2d Landau-Ginzburg 
theories and from the point of view of continuum theory they should be identified with 
the twisted Landau-Ginzburg topological matter interacting with topological gravity [5 

(Of course, the basic vectors of the pure (p, 1) model corresponding to monomial 
potential V p can be obtained by setting t\ = ... = t p = 0, t p+ i = -^j). As a solution to 
string equation this deformed case differs only in analytic continuation along first p times. 

These topologically deformed (p, 1) models as we already said preserve all topological 
properties of (p, 1) models. Indeed, according to shifting of first times ti,...,t p+ i is 
certainly not enough to get higher critical points. To do this one has to obtain t p+q ^ 0, 
but this cannot be done using above formulas naively, because it is easily seen from 
definition ( |223| ) of p-times, that t k = for k > p + 2. 



9 General scheme and ^-duality 



The above scheme has a natural quasiclassical interpretation. Indeed, the solution to 
(jp, 1) theories given by the partition function (f76| ) can be considered as a "path integral" 



representation of the solution to Douglas equations 

[P,Q] = 1 (240) 

where P and Q are certain differential operators (of order p and q) respectively and 
obviously p — th order of P dictates p-reduction, while q stands for q — th critical point. 
Quasiclassically, (|240|) turns into Poisson brackets relation [|36|, |37] 



{P,Q} = 1 (241) 

where P(x) and Q(x) are now certain (polynomial) functions. It is easily seen that 
the above case corresponds to the first order polynomial Q(x) = x and the p-th order 
polynomial P(x) should be identified with W(x) = V'(x). Thus, the exponentials in (|76|), 
and (|238|) acquire an obvious sense of action functionals 



S pA (x, ii) = -V(x) + xWfji) =- dy W(y)Q'{y) + Q(x)W(n) 

Jo 

V 

W(x) = V\x) = x p + J2 v kX 



fe-i 

k=l 

Q(x) = x u (242) 



and we claim that the generalization to arbitrary (p, q) case must be 

S W)Q = ~ f dy W(y)Q'(y) + Q[x)W{ji) 

•J 



W(x) = V'(x) = x p + ^2 v kX 



fe-i 

UkX 

fc=l 

9 



Q( x ) = x q + J2v k x k - 1 (243) 



k=l 



Now the "true" co-ordinate is Q, therefore the extreme condition of action (|243|) is still 



W{x) = W{ii) (244) 



14 



For example: W(z) = z 2 + t x , Q(z) = z, then for ( |24lD 



having; solution, and for extreme value of the action one gets 



S w ,q\ x= , = I dy W'(y)Q(y) 



o 



where jjP = W(fi) and 



E (245) 

fc=— oo 



t k = tf ' Q) = -—P—Res W^dQ . (246) 
k(p — k) 



We should stress that the extreme value of the action (|243|) , represented in the form 



( |245| ), determines the quasiclassical (or dispersionless) limit of the p-reduced KP hierarchy 
p]| , P7[ | with p + q — 1 independent flows. We have seen that in the case of topologically 
deformed (p, 1) models the quasiclassical hierarchy is exact in the strict sense: topological 
solutions satisfy the full KP equations and the first basic vector is just the Baker- Akhiezer 
function of our model ( [75] ) restricted to the small phase space. Unfortunately, this is not 
the case for the general (p, q) models: now the quasiclassics is not exact and in order 
to find the basic vectors in the explicit form one should solve the original problem and 
find the exact solutions of the full KP hierarchy along first p + q — 1 flows. Nevertheless, 
the presence of the "quasiclassical component" in the whole integrable structure of the 
given models can give, in principle, some useful information, for example, we can make a 
conjecture that the coefficients of the basic vectors are determined by the derivatives of 
the corresponding quasiclassical r-function. 

Returning to eq. ( |246|) we immediately see, that now only for k > p + q + 1 p-times 
are identically zero, while 

f = f (w,Q) _ P (247] 
i p+q -i p+q + y V^t) 

and we should get a correct critical point adjusting all {t k } with k < p + q to be zero. 
The exact formula for the Grassmannian basis vectors in general case acquires the form 

<M/i) = [W'Gu)] 1/2 exp(- S w , Q \ x=t ) J dM Q (x)f t (x)e* V S w , Q (x,n) (248) 

where <1Mq(x) is the integration measure. We are going to explain, that the integration 
measure for generic theory determined by two arbitrary polynomials W and Q has the 
form 



by checking the string equation. For the choice (|249 ) to insure the correct asymptotics 



of basis vectors 4>i{fj) we have to take ft(x) being functions (not necessarily polynomials) 
with the asymptotics 

/^-^(l + CKl/x)) (250) 

To satisfy the string equation, one has to fulfill two requirements: the reduction con- 
dition 

W{riUri = Y.C l ^) (251) 

j 

and the Kac-Schwarz (|118|) operator action 



A^V^) = E^<Mm) (252) 

with 

A (W, Q) ^ ^WO) M _^|. [JV WO) (Al)] -i = 
1 d lW"(/x) . 



W'(n)dn 2W'(^) 2 
iV™^) ee [W'i^eM- Sw, Q \ x= ,) (253) 

These two requirements are enough to prove string equation. The structure of action 
immediately gives us that 

A^Uri = N^ Q \ri J dM Q (z)Q(z)f l (z)exp S w>Q (z,fx) (254) 



and the condition (|252|) can be reformulated as a Q-reduction property of basis {fi(z)} 

Q{z)f i {z) = Y,^fi{ z ) ( 255 ) 

Let us check now the reduction condition. Multiplying 0i(/i) by W(fjj) and integrating 
by parts we obtain 

= N^\n) J dM Q (z)ft(z)^^[exp Q(z)W(n))exp[- J* dy W{y)Q'{y)] = 
= -NWMfr) [ dM Q {z) bMSw,q{z^)\ ( t^#- - - W{z)\ f { (z) = 



Q'{z)dz 2Q'{zf 



Therefore, in the "dual" basis {fi(z)} the condition (31) turns to be 



A«^/,(z) = -£<Vi(*) (257) 
with A( Q ' w \j£ A {w ^) being the "dual" Kac-Schwarz operator 

- 4,0 '"' , -Q^l-l^W- ™ < 258 > 

The representation ( 248Q , (|249 ) is an exact integral formula for basis vectors solving 
the (p,q) string model. It has manifest property of p — q duality (in general W — Q), 
turning the (p, g)-string equation into the equivalent (q, p)-string equation. 



Now let us transform ( 248 ), ( 249 ) into a little bit more explicit p — q form. As before 
for (p, 1) models we have to make substitutions, leading to equivalent KP solutions: 



/F = W(n), z q = Q(z) 



Then we can rewrite (248) as 



Up) = [p/F^pexp (-E^J / d~z[q~z^ l } l l 2 h{~z)e W S w<Q (z,fl) 
where action is given now by 



S WiQ (z,jl) = - \ T dyqy^Wiviy)) 
Jo 



+ z q p? 



P+Q _ 

fc=i 

In new coordinates the reduction conditions are 



(259) 



(260) 



(261) 



and for the Kac-Schwarz operators one gets conventional formulas ]2"3|, |2~4] , |TB 



i ft -n - 1 1 1 p+q 



ppP 1 dp, 2p pP p£r{ 
1 d q - 1 1 1^ 



1.4- zk-q 



(262) 



where for (q,p) models we have introduced the "dual" times: 



(Q,W) 



Res Q 



(264) 



k(q — k) 



in particularly, t p+q = —H p+q = — ^- . Now string equations give correspondingly 



By these formulas we get a manifestation of p — q duality if solutions to 2d gravity. 
As a main result of formulas presented above, one may conclude that a generic solution 
to c < 1 2d gravity should be described by two functions W(x) and Q(x). 



Now we will make some remarks why this model could be considered as an attempt of 
constructing a string field theory or effective theory of string models. It is necessary to 
point out from the beginning that by string field theory we would mean more than a con- 
ventional definition as a field theory of functionals defined on string loops - it must rather 
mean a sort of effective theory which gives all the solutions to classical string equations 
of motion (2d conformal field theories coupled to 2d gravity) as its vacua and allows us to 
consider all of them on equal footing (within the same Lagrangian framework) and maybe 
even describe the flows betveen different string vacua. Of course, it should reproduce 
perturbation expansion around any of these vacua. In this sense, the conventional string 
field theory was not true effective model, because it contained a fixed set of variables 
which correspond to a concrete vacuum (say, 26 free scalar fields). Therefore, it doesn't 
have even a priori a possibility to make a flow to another classical solution (maybe only 
except for some simple change of a background), i.e. conventional string field theory 
might describe only some small perturbation around given classical solution in terms of 
the coordinates equivalent to the matter variables in the Polyakov path integral ([!]). 

Moreover, in conventional approach even the perturbative expansion was ill-defined 
due to the presence of tachyon in the spectrum or in other words due to the instability of 




(265) 



10 String field theory and c — » 1 limit 



make sense to the non-perturbative description, presented above, and existing up to now 
only in the case of non-critical and moreover "non-tachyonic" strings. That means that 
what we know up to now is the only case of "highly-noncritical" string models where the 
total matter-gravity central charge 

^matter ^gravity 26 (266) 

is "dominated" by the contribution of 2d gravity. This is far from the case of critical string 
(cmatter = 26) and close to the case of pure two-dimensional gravity. Unfortunately, up 
two now this is the only region where it is possible to formulate string theory consistently 
and at least put the question what is the internal principle which might allow one to 
choose dynamically a string vacuum. 

Now, the above scheme is a string field theory in a sense it is highly related to the 



underlying module space |33], . That is actually the main difference with ordinary field 



theory: the geometrical sense of construction above was motivated by cell-decomposition 



of module space [§g, |47|, |33], |34 



So, at the moment we have a sort of theory, describing various (p, q) models coupled 
to 2d gravity beyond the perturbation theory, and technically advocating to the equations 
in the space of coupling constants, which determine the generation function for all the 
correlators. Such function (being the logariphm of the r-function cannot be defined 
globally in this space, and around each "critical point" it has some sort of perturbative 
definition, which should reproduce the original "first- quantized" theory (U). However, 
moving in this space from one solution to another naively one will meet the divergences 
of such expansion, then it is necessary to continue analytically, and the piece we would 
get in such a way is exactly a nonperturbative correction. For simplest (p, 1) topological 
theories, this scheme can be described at the language of effective ("path") or rather 
matrix integral with almost trivial integration measure, the exact integral formulas for 
generic case involve more complicated structures. 

This scheme, in principle, should also be true for the "barrier" c = 1 solutions. How- 
ever, a generic c = 1 "phase" is rather complicated theory, and naive one-matrix formu- 
lation leads only to some highly restricted c = 1 cases. 



the determinant form of Penner model ||6, |47], [18] partition function implies already that 
for fixed values of times it is a Toda lattice tau-function in the sense of and allows us to 
apply to this case the the Toda theory representation for Generalized Kontsevich models. 
Indeed, the solution to the Penner model 



Z ~ det U 



(a) 



with 



Hff = T(a + t+j-l) 



(267) 



(268) 



is nothing but a specific case of GKM. 



Now one can easily introduce positive- and negative-times dependence in (|267|) and 
then reconstruct $^ (z) from (162) Indeed, 



=^ = T{a-l+i+j) 



immediately gives 



dy 
o y 



e ~VyOl-l+i+j 



dy 
o y 



(269) 



(270) 



which is a sort of GKM-like representation. The difference with more common situation 
for c < 1 is in the definition of the contour in (|27C|) and also in the fact that z-dependence 
is trivial, because integral is easily taken with the result 

V(a + i) 



and 



9 1 A<*), 



(z-1) 



oA a ) , 



a+i 



Introducing negative times, one gets P0 



■ Tja + i+j) 



(271) 



(272) 



z|T_ p ) = z' 



exp -£T_ 



p>o l k 

where Pk[T p ] are Schur polynomials (exp J2T p z p = z n P n [T p }). or simply 
Z_i ~ f DY exp TrZY + aTrloaY + V T,TrY- k 



(273) 



f274) 



with 



T +k = yTrZ k (275) 
k 



and amazingly this formula turns to be consistent with the calculations of the tachyonic 
amplitudes |45[ . 

Let us now make more comments on c = 1 situation. From basic point of view we need 
in generic situation to get the most general (unreduced) KP or Toda-lattice tau-function 
satisfying some (unreduced) string equation. In a sense this is not a limiting case for 
c < 1 situation but rather a sort of "direct sum" for all (p, q) models. This reflects that in 
conformal theory coupled to 2d gravity there is, in a sense, less difference between c < 1 
and c = 1 situations than this coupling. 

However, there are several particular cases when one can construct a sort of direct 
c — > 1 limit and which should correspond to certain highly "degenerate" c = 1 theories. 
From the general point of view presented above these are nothing but very specific cases of 
(p, q) string equations, and they could correspond only to a certain very reduced subsector 
of c = 1 theory. 

Indeed, it is easy to see, that for two special cases p = ±q the equations ( [25 1| ), ([25 2|) 
can be simplified drastically, actually giving rise to a single equation instead of a system 
of them. Of course, these two cases don't correspond to minimal series where one needs 
(p, q) being coprime numbers. However, we still can fulfill both reduction and Kac-Schwarz 
condition and these solutions to our equations using naively the formula for the central 
charge, one might identify with c = 1 for p = q and c = 25 for p = —q. 

Now, the simplest theories should be again with 5 = 1. For such case "c = 1" turns to 
be equivalent to a discrete matrix model |25| while "c = 25" is exactly what one would 
expect from generalization of the Penner approach |50 , [55 1. Indeed, taking in general 
non-polynomial functions, like 

W(x) = x~ p 

Q(x) = x p (276) 
the action would acquire a logariphmic term 

S- Pt0 = -(3logx + (277) 



while equations (|251|) , ( |252| ) give rise just to rational solutions. It is very easy to see that 
(3 = 1 immediately gives the Penner model in the external field, which rather corresponds 
to "dual" to c = 1 situation with matter central charge being c maiter = 25 with a highly 
non-unitary realization of conformal matter 0. 

On the other hand, p = q = 1 solution is nothing but a trivial theory, which however 
becomes a nontrivial discrete matrix model for unfrozen zero-time. Moreover, these par- 
ticilar p = ±q solutions become nontrivial only if one considers the Toda-lattice picture 
with negative times being involved into dynamics of the effective theory. On the contrary, 
we know that c < 1 (p, g)-solutions in a sense trivially depends on negative times with the 
last ones playing the role of symmetry of string equation |23J . It means, that we don't yet 
understand enough the role of zero and negative times in the Toda-lattice formulation. 



11 Conclusion 

Let us, finally make several concluding remarks. We tried to demonstrate above that using 
additional structure of integrability arised in the "discretized" formulation of 2d gravity 
one might go far from what is known in string perturbation expansion. Namely, a generic 
(p, q) solution to the c < 1 case can be described as a solution to the integrable system and 
its explicit form can be studied using various integral representations for these "stringy" 
solutions. In the most simplest case of topological theories, these integral formulas can be 
combined into the matrix integral giving rise to a matrix "zero-dimensional" field theory 
description. 

Much more complicated is the case of d > 1 which at the moment seems to be almost 
unclear from the point of view be discussed below. One might only hope that there exists 
a possible generalization concerning non-Cartan flows and formed by them an integrable 
system. One might also hope that such integrable systems are related to the matrix-like 
integrals with a nontrivial integration over "angle variables", i.e. when in matrix integral 
(HD the action is no longer the trace of function of a matrix, but is rather complicated 



objects. Such type of systems were discussed in the literature [53] and the technique 



15 This c = 1 - c = 25 duality might be also connected with the known fact that there exists a Legendre 
transform between the Gross-Klebanov solution to c=l matrix model and the Penner model 



was mostly based on the more complicated operations with the Itzykson-Zuber integrals 
[plfl . This might be connected to a sort of "group-theoretical" r-functions, to be discussed 
elsewhere [|54] . 
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